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Abstract. Based on our previous analysis Q of the anisotropic integrable chain 
consisting of spins s = | and s = 1 we compare the dispersion relations for the sectors 
with infinite Fermi zones. Further we calculate the speeds of sound for regions close to 
' sector borders, where the Fermi radii either vanish or diverge, and compare the results. 
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1. Introduction 



The integrable spin chain XXZ(~, 1) constructed in 1992 by de Vega and Woynarovich 

shows a rich physical structure with different ground states depending on the 
anisotropy parameter and the two coupling constants. In our previous papers pj, 0, 

01 and || we have analyzed this structure and calculated several important quantities. 
In this paper we deal mainly with dispersion relations and speeds of sound. 

The paper is organized as follows. After having reviewed the definitions in section 

2 we found it necessary to remind the reader results of our last paper [|J concerning the 
phase diagram. In the same section the dispersion relations for all sectors with infinite 
Fermi zones only are compared. In section 4 we present the calculations for the speeds of 
sound obtained either by power expansion or Wiener-Hopf technique for small or large 
Fermi radii respectively. 

Our conclusions are contained in section 5. Some useful definitions are compiled in 
an appendix. 



2. Description of the model 

We refer the reader to papers |2| and || for the basics of the model. 
Our Hamiltonian of a spin chain of length 2N is given by 

H( 1 ) = cH( 1 ) + cH( 1 ), (2.1) 

with the two couplings c and c. The anisotropy parameter 7 is limited to0<7<7r/2. 
For convenience we repeat the Bethe ansatz equations (BAE), the magnon energies and 
momenta and the spin projection. 

/ sinh(A j + q) smh(Xj + ij) \ N _ |4 sinh(A, - A fc + 27) ^ M 

ysinh(Aj — i%) sinh(Aj — ry) J M sinh(Aj — A^ — 27) ' 

E = cE + cE, 
- = _« 2sin 7 
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-J cosh 2\j — cos 7 
• u 2 sin 27 



E = -Y — 5— r 1 , (2.3) 

cosh 2Xj — cos 27 
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S z = — -M. (2.5) 



M 



'sinh(Aj + ( sinh(Aj + ry)' 



Dispersion relations and speeds of sound 



3 



o 


(1 +)°° (2 +)°° 
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(1 )oo (1 + )/,0 


II 


(1,-)°°, (1,+)°° 


III 


(1,-)°°, (1,+)*°° 


IV 


(I,-) 00 


V 


(i,-r, (2,+p 


VI 


(1,-)/'°°, (2,+)°° 



Table 1. All sectors appearing in the phase diagram. Upper indices indicate infinite 
and finite Fermi zones. In the latter case, the second index distinguishes, wether the 
filling starts at A = or A = oo. 



3. Phase structure and dispersion relations 

In section 4 of [J] we have carried out the complete analysis of the ground state structure 
based on the solution of the thermodynamic Bethe ansatz. We will now review the main 
results. 

Using different sets of special anisotropy points we were able to prove that (for 
< 7 < 7r/2) only three kinds of strings occur in the ground state. The picture changes 
remarkably at 7 = 7r/3 and 7 = 2ir/5. All possible sectors are compiled in table 1 taken 
from To keep our notation consistent we found it useful to call the sector already 
discussed in sector 0. In the (c, c)-plane the sector borders are straight lines starting 
from the origin. Moving counterclockwise from the c-axis the sectors follow in sequence 
0, I, IV, V if < 7 < tt/3. For tt/3 < 7 < 2vr/5 the sequence is 0, II, I, IV, VI and for 
27r/5 < 7 < tt/2 one has 0, II, III, IV, VI. In case 7 = 2tt/5 sector II is followed by IV 
and neither I or III occur. For details see figure 1 in 0. 

To compare the dispersion relations in the sectors 0, II and IV it is convenient to 
work with the functions 

sinh(7rA/ (3) 

P{\a,P) = arctan — — (3.1) 

cos(7ra/2p) 

and 

4tt cos(7ra/2/3)cosh(7rA//5) 

g{A,a,p) — — - j—r. [6.Z 

(3 cosh(27rA/p) + cos(7ra/2p) 

While p is odd, g is an even function of A. Both are connected by 

-^p(A, a, (3) = ~#(A, a, (3). (3.3) 
By straightforward calculation one may prove the relation 



7T 



sin I — +p(\,a,(3) 



, sin 2 p(A,a,/5) 
cos 2 p{\,a, (3) + — r = — g(\,a,/3). (3.4) 



.2 J \ cos A { r Ka/2[3) 2n 

We remark that in || we used a function g(X, a) related by 

^(A,a,7r-7) =g(X,a). (3.5) 
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Let us first consider sector 0. The results are given in |2[ taking into account our 
definition of momentum in equation ( |2.4| ). One has two elementary physical excitations, 
the holes in the distributions of the (1, +)- and (2, +)-strings respectively. We label them 
by index i, i = 1, 2. 

PUA = A1) = | + ^(A,0, 7 ), (3.6) 

4(A) = <*\g{\ 0, 7 ) (3.7) 

with c 1 = c and c 2 = c. 

From equation (|3.4|) the dispersion relations follow 



4«) = -c^. (3.8) 
7 2 

One may say that the dispersion relations diagonalize in the coupling constants. This 
fact changes in the other two sectors. 

Next we look at sector IV which contains the case of negative couplings. The 
exact border lines are given in ||. The elementary physical excitations are holes in the 
distribution of the (1, — )-strings. We quote the result from our paper |J, equations 
(4.5) and (4.6) 

PhW = Y + \p{\ 0, vr - 7 ) + p(A, 7, vr - 7), (3.9) 

= ~9{\ 0, 7T - 7) - cg{\, 7) tt - 7 ). (3.10) 

The dispersion relation cannot be obtained in a closed form. One has to invert the 
monotone function (|3.9|) and to substitute it in equation (|3.10|) . Numerical calculations 
show an sine-like behaviour, see the broken lines in figure 2 of 0]. 

In [|TJ we found the new sector II, existing only for 7 > tt/3 (for the exact border 
lines see next section). Here the excitations are (1,+)- and (1, — )-string holes. The 
dispersion law for the (1, — ) -string holes is simple: 

pjt-(A) = £ + ^(A,0 l 7r-2 7 ), (3.11) 
4"(A) = -^(A,0,7r-2 7 ), (3.12) 

= - 'I Sil ?^ (3-13) 

7T — Z7 L 

To write down the case of (1, +)-strings it is necessary to introduce a new function di(X) 

A ( X ) = 1- r r ^A cosh(^(7r - 3 7 )/2)^ 
1K> 27t7-oo 2cosh(cu(7r-2 7 )/2)cosh(^ 7 /2) K ' 

already used in Jl|. We further define another function ei(A) by 

A 6l (A) = d x {\) and ei(0) = 0. (3.15) 
uA 




Figure 1. Dispersion relations for the holes in the (1, +)-string distributions in sector 
II for c = —5c = 1 and different 7. 



^ + = ^ + ^(A,0, 7 ) + 27r ei (A), (3.16) 

4 + = c^(A,0, 7 ) + 27Tgd 1 (A). (3.17) 

We have calculated the dispersion law numerically and show the results in figure [T]. It 
is interesting to look closer at the point 7 = 2ir/5, where the function di(X) simplifies 
and becomes l/47rg(A, 0,7), resulting in the simple dispersion law 

1+ 27TC + C . 1+ 

e h + = — smp h + . (3.18) 

7 2 

All dispersion laws possess a symmetry connected with the properties of the functions 
p h (X) and e h (X): 

PhW=Po+P(X) and p(-\) = -p(\), (3.19) 

e h (X)=e h (-X). (3.20) 
One then derives easily 

£h(Ph) = £h(2po ~ Ph)- (3.21) 
In addition one has < ph < 2p . The constants 2p are different, so in sector they 
are 7r/2 for both excitations, the same is true for the negative parity excitation in II. 
The other (positiv parity) excitation in II has 2p = it, while in IV one has 2p = 3ir/2. 
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There is a general rule for the sum of those constants in every sector being 37r/2 (due 
to antiferromagnetic ground state), but excitations containing more than one magnon, 
e.g. (2, +)-strings, must be taken with their appropriate multiplicity. 

4. The speed of sound in different sectors 

In this section we shall present some calculations for the speed of sound in sectors with 
one finite Fermi radius (the other one stays infinite) where this radius is either small or 
large, so that an approximation can be made. On the other hand, that gives just the 
values of the speed near to the border lines of the above mentioned sector which might 
be of interest for understanding the physics of the model. We shall limit our calculations 
to the case c > 0, c < 0. Thus, the sectors in question are I and III. Sectors V and VI 
could be treated in an analogous way. 

We start with sector I. The equations of the thermodynamic Bethe ansatz can be 
taken from |] (equation (4.5)). 

e 1+ (A) = -2ttcsi(A) - 2ncd l (X) + d x * e+ + (A), (4.1) 

ei_(A) = 27rg S2 (A)-s 2 *e+ + (A). (4.2) 

Here Ei{X) are the dressed energies. Their lower index specifies the excitation, while the 
upper stands for the positive and negative parts respectively. The convolution is defined 
by 

POO 

a * 6(A) = / d[ia(\ — fi)b(fi). (4.3) 



(For Si(A) and s 2 (A) see the appendix.) 

Restricting 7 by < 7 < 2ir/5 we consider the case 

0</3 = c--^<l (4.4) 
2 cos 7 

with the usual definition 7 = 717/2(7? — 7). 

From equation (3.13) from Q we know that (3 — describes the border line between 
sectors IV and I. 

To employ the method of power expansions in the (small) Fermi radius b of the 
(1, +)-strings, it is useful to reformulate equations ([O]) and ([12]) in the way 



£ i+(A) = -c#(A, 7T - 27) - cg(X, 7T - 37) - ^-g(X, 7T - 37) * e 1+ (A), (4.5) 

ei_(A) = 2ttc Si (A) + cg{X,j) + -^(A, 7) * Cf(A). (4.6) 

Equation Q4.5Q determines the Fermi radius b, in lowest order the convolution term be 
neglected and b is given by = 0. The derivative de/dp is calculated the way (see 

equation (5.1) in [|TJ) 

de de dX de 2 



dp dXdp dXe(X) 



(4.7) 

c=c=l 
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The speed of sound is then given if A is put equal to b. 

After rather long technical calculations we have obtained 



47T 1 /2cos7 - 1 , ml/2 f A 

«i+ = — Vo ir—rWr- 4 

7r — 7 sm 7 y 2 cos 7 + 1 



Equation gives the same way 



2^ /ri ~ ,\ 8 . _ / 2 cos 7 — 1 , ..I/O , . 

vi- = (2cos 7 -l)c+ sm T W- ^—(c/3) 1 / 2 . (4.9) 

7r — 7 7r — 7 y 2 cos 7 + 1 

Here the first term is just the limiting value on the border of sector IV. 

Next we consider 7r/3<7<27r/5 and look for the border of sector I with sector II 
(see [Q, equation (4.10)). It is given by 

< a = |c| tan f^J - c < 1. (4.10) 

Then equation (|4.1| ) is of Wiener-Hopf type and can be solved explicitely for large b 
giving 

2tt 37 — 7r a ,. ... 

7 vr-271+tan(g) 

From equation ( |4.2| ) one has to extract the asymptotics of £i-(A) for large A (Fermi 
radius is infinite). The first term on the RHS gives the limiting value of on the 
border of sector II (it is non- vanishing) . Remarkably for a > the asymptotics is 
always given by the second term with the convolution. So the result is of the same order 
in a as V\ + : 

2n 37 — 7r a , d . 

vi- = 7~2\- 4 - 12 

7 7 l+tan(g) 

Therefore v\- shows a discontinuity on the border I /II. 

After sector I we consider sector III. Then one has 2ir/5 < 7 < ir/2. We start again 
with the border to sector IV. 

< ft = c-2cos7C< 1. (4.13) 

The Fermi radius is determined again from equation ( f4.5|) , but now instead of power 
expansion we have to solve a Wiener-Hopf problem due to the fact that b is large. The 
calculation is straightforward after realizing that for the positive frequency part of the 
inhomogenous function one has to take into account the dominant two poles (instead 
of one). As already mentioned in the fact that filling of the (1, +)-strings starts at 
infinity causes the appearance of two velocities for them taken at A = b and A = 00. 
While the calculations for both are different the results coincide in first order. 

The same happens for calculated from the asymptotics of equation ( |4.6| ). We 
remind the reader that v\- vanishes at the border. Finally we have obtained 



m (2) 27T & 
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We have no explanation for that striking fact, that all velocities are equal in first order. 
Nevertheless we are sure that the effect does not survive in higher order. 

To close that part we have to consider the border line with sector II (see Jlj] , equation 
(4.8)). 

< a' = |c|2 7 rfi(0) -c< 1. (4.15) 

The Fermi radius b is determined by equation ( [4.1|) using power expansion. The 
calculation is straightforward, but rather tedious because of the presence of the function 
g?i(A). But the order of magnitude can be estimated easily. The two velocities ?4_ and 

(2) 

v\_l (taken at A = oo) obtain corrections to their (constant) limiting values from sector 
II of the order \fd . The third velocity v j+ (taken at A = b) is of order \[a! '. 

We did not consider the border of sectors I and 0for0<7<7r/3 because both 
sectors are separated form each other by a singular line with a highly degenerate ground 
state which is different from the other lines considered above. 

5. Conclusions 

First we make some general comments on the character of the behaviour of the velocities 
considered. We parametrize the coupling constants by two parameters, one c, setting 
the energy scale (below put to one), the other one an angle (p, which describes their 
ratio. 

c = c sin (p, 

c = ccosip. (5.1) 

In general from our results in section 4 one can see that the behaviour near the sector 
borders is either linear or a square root in cosine/sine- functions of (p. The fact which 
one is chosen is determined by the radius b. Large b produces Wiener-Hopf equations 
and linear behaviour, small b causes power expansions and a square root behaviour. 
Considering sector I only, one could superficially conclude that this depends only on 
the sector on the other side of the border. Our analysis in sector III shows that the 
situation reverses with respect to sector I. 

In just one case we have found a discontinuity of one of the velocities when a sector 
border is passed (see equation ( |4.12| )). It is the case for t>i_ on the border I/II. Nothing 
of this kind happens for sector III, because for small a' we have power expansion, non- 
capable of producing such an effect. On the other hand (for small f3') V\_ vanishes on 
the other side of the border in sector IV. In general, one would expect a discontinuity 
only when a singular line is crossed which means that an excitation with infiniteFermi 
zone disappears at once (see below). It is rather instructive to compile all analytically 
calculated velocities in one picture at the special point 7 = 2tt/5 which we have done 
in figure |2|. 

For convenience, we have given the analytical shape in the various sectors in table 
H Looking at figure || one clearly sees two effects. First the two singularity lines at 
tp = 7r/2 (the vanishing of the (2, +)-strings) and at p> = 3ir/4 (the vanishing of the 



Dispersion relations and speeds of sound 



9 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 2. The speeds of sound as functions of the angle tp for 7 = 2ir/5 and those 
sectors where all speeds of sound are known. 



sectors 





II 


IV 


V 


borders 


0, tt/2 


7r/2,37r/4 


37r/4, 27r — arctan4 


27r — arctan 4, 27r 


«i+ 


5 simp 


5/2 sin ip + 5/2 cos yo 








V2+ 


5 cosy? 










V\- 





— 10 cosy? 


—5/3 sin ip — 5/3 cos ip 





Table 2. The speeds of sound and the sector borders as functions of the angle ip for 
7 = 2tt/5. 



(1, +)-strings) cause discontinuities for the speed of the remaining excitation. While the 
first line is universal for arbitrary 7, the second one is peculiar to 7 = 2ir/5. For other 
7 that line is replaced by sectors I or III. 

Second one sees the two conformal points where two velocities intersect, one at 
ip = 7r/4 and the other at <p = it — arctan 5 |4]]. 

Our calculations in the previous section allow now to look closer at what happens 
for 7r/3 < 7 < 27r/5 when the above mentioned singularity line is replaced by sector I. As 
far as the limits of sectors II and IV are concerned the only change is that t>i_ no longer 
vanishes at the border of sector IV. Now we make use of equations (|4.9|), ( |4.10| ),( |4TTl"| ) 
and ( |4.12| ). They tell us first that both velocities must reach their maximum (at different 
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<p) within sector I and second that there is (at least one) point of intersection in the same 
sector. In that point the finite size corrections must exhibit a rather simple structure 
like in §. 

We conjecture that there is no such effect in sector III but a strict prove requires 
further calculations. In addition, our calculations show another interesting effect, at 
(3' — all velocities vanish. 

Finally we remark that we expect the theory at 7 = 2ir/5 and both velocities equal 
to be highly symmetric and of particular interest. Up to now it is only known that its 
central charge is 2. We were not able to calculate the operator dimensions analytically, 
because contrary to what happened in case of paper || it requires the knowledge of the 
dressed charge || and hence the explicit factorization of the kernel. 



Acknowledgment 

The authors would like to thank H. Frahm for helpful discussions. 



Appendix 



Sl(A) 2 7 cosh(vrA/7 ) 
s 2 (A) ' 



2(tt - 2 7 ) cosh(7rA/(7r - 2 7 )) 

f\ \ f° «,a coshwa/2 

g(X,a)= e — — — — du 

J -00 coshu;(7r — 7J/2 
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Figure and table captions 

Figure 1. Dispersion relations for the holes in the (1, +)-string distributions in sector 
II for c = —5c = 1 and different 7. 



Figure 2. The speeds of sound as functions of the angle ip for 7 = 2n/5. 



Table 1. All sectors appearing in the phase diagram. Upper indices indicate infinite 
and finite Fermi zones. In the latter case, the second index distinguishes, wether the 
filling starts at A = or A = 00. 



Table 2. The speeds of sound and the sector borders as functions of the angle ip for 
7 = 2vr/5. 



